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In this article we consider the Hamiltonian dynamics of systems of fermions 
and derive the time-dependent Hartree-Fock equation in the mean field limit. 
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We follow the approach of Spohn, who derived a mean field dynamical equa- 
tion (the time-dependent Hartree equation) for mean field systems of dis- 
tinguishable particles, remarking at the time that "the convergence of the 
mean field limit with statistics included is an open problem" . 

In Spohn's theory the initial A'^-body density operator D^^ is assumed 
to be a product state D^^ , i.e., the particles are statistically indepen- 
dent and identically distributed. The mean field limit is investigated in 
the Schrodinger picture, where D]\j{t) obeys the von Neumann equation 

ihjD^it) = [LvDNit)] + ^ E i^^J^Dr^it)] (1) 

l<j<N ^<i<j<N 

with Vij denoting the two-body potential V acting between the i*^ and j*'^ 
particles and [ , ] denoting the commutator. The limit as N — > oo of the n- 
body density operator D^-nit) is shown to converge to D{t)®^ , where D{t) 
obeys a time-dependent Hartree equation. (The subscript appearing in 
D^-.n is our notation for the n*'* partial trace, defined in equation (Q) below.) 
Spohn's ideas have been generalized to open systems in Q. There are other 
theories of quantum mean field dynamics, e.g., the algebraic theory of 
but to our knowledge the problem of including quantum statistical effects 
remains unsolved. 

The problem is that Fermi-Dirac or Bose-Einstein statistics constrain the 
possible initial condition of (||) to have the appropriate symmetry, which is 
typically inconsistent with the product form D®^ . An A'^-body density 
operator with Fermi-Dirac symmetry can never have the form D®^ and a 
Bose-Einstein density operator can only have the form D®^ if D is a pure 
state (i.e., if the system of bosons is in a condensed state). The remedy for 
this problem, for fermions, is to replace the hypothesis that the initial state 
be a product state with a hypothesis that is consistent with Fermi-Dirac 
statistics, e.g., that the initial states are Slater determinants. 

The role of the factorization hypothesis D]\f{0) = D®^ is to permit the 
closure of the BBGKY hierarchy by setting the two-body state Dn-2 equal to 
D®D. Closing the hierarchy this way results in the time-dependent Hartree 
equation. This kind of closure hypothesis is implicit in the Stosszahlansatz 
that leads to Boltzmann's kinetic equation for gases [Q]. Kac noted that, 
for Boltzmann's equation, the factorization /7V:2 = / <^ / is only realized 
in the limit — > oo, and he called this behavior the Boltzmann property 



[12, 13|. Later authors [14, 11, 17[ developed Kac's ideas; what is now called 
the propagation of chaos is an important tool in rigorous kinetic theory [^, 
15, 0]. We have noted that Boltzmann's closure Ansatz is inconsistent with 



2 



the Pauli Exclusion Principle, and needs to be replaced by another closure 
Ansatz when the particles are fermions. The novelty of our approach consists 
in replacing the condition of asymptotic independence of the particles by 
a condition that describes the correlations of Slater determinants. This 



condition, called Slater closure is defined in Definition 2.1 below. 

Assuming that {Di\f{0)} is a sequence of initial states for (|l]) that has 
Slater closure, we can prove that {Diy{t)} has Slater closure for all t > 0. 
This phenomenon could be called the propagation of Slater closure because 
it is like the "propagation of chaos" mentioned above. Since {D]y{t)} has 
Slater closure, the two-body density operator Djsi:2{t) is approximately equal 
to {Di^i-i{t)(^DN-i{t))T,2 when N is large, where S2 is the two-body operator 
defined by 

'^2{x(^y) = x®y — y®x. 

Substituting {DN:i{t) ^ DN:i(t))T.2 for DN:2{t) in the BBGKY hierarchy 
leads one to conjecture that, when is large, the single-body density oper- 
ator should nearly obey the time-dependent Hartree-Fock (TDHF) equation 

ih^F{t) = [L, F{t)] + [V, (F(t)®F(t))S2].i 
dt 

F(0) = Dn:i{0). 



Theorem confirms this conjecture. 

Theorem states that the distance in the trace norm between DN:i{t) 
and the corresponding solution F{t) of the TDHF equation tends to as A 
tends to infinity. The trace norms of D]\}-i(t) and F{t) are separately equal 
to 1, so it is significant that their difference Diy-i{t) — F{t) converges to 
in the trace norm. A crucial detail of the proof is Lemma |5.1| , which states 
that the operator norm of Djy-i tends to if {-Dat} has Slater closure. Much 
of the rest of the proof lies in bounding the trace norm of Diy-i(t) — F{t) by 
an expression involving the operator norm of D]\f:i{0). 

The use of the trace norm to measure the distance between two density 
operators is quite natural. A density operator D corresponds to a quantum 
state through the assignment B 1— > Tr{DB) of expectation values to bounded 
observables B. Thus, two density operators D and D' are within e of one 
another in the trace norm if and only if they correspond to quantum states 
that give expectations differing by no more than e for all observables B with 
ll^ll < 1- 

In this article, we assume that the two-body potential F is a bounded 
operator. Recent work on the time-dependent Schrodinger-Poisson equation 
[0, suggests that it may be possible to prove a theorem similar to our 
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Theorem |3.1| when V is the Coulomb potential. This work shall be published 
in a separate paper. 

The rest of this article is organized as follows: the next section dis- 
cusses fermionic density operators and defines Slater closure. The N-particle 
Hamiltonian and the associated time-dependendent Hartree-Fock equation 
are described in Section 3. That section concludes with the statement of 
our main result, Theorem 3.1, whose proof spans Sections 4, 5, and 6. Sec- 
tion 7 is an appendix relating the operator form of the TDHF equation, 
used throughout this paper, to the formulation of the TDHF equation as 
a coupled system of wave equations, which may be more familiar to some 
readers. 



2 Fermionic density operators and Slater closure 

Let M he a Hilbert space, supposed to be the space of wavefunctions for a 
certain type of quantum system (a "component" or "particle"). Then the 
Hilbert space of wavefunctions for a system consisting of N distinguishable 
components or particles of that type is Mn = M"^^ . If the components are 
not distinguishable, but obey Fermi-Dirac statistics, then the appropriate 
Hilbert space of wavefunctions is the antisymmetric subspace An ^n- To 
define this subspace, it is convenient first to define unitary transposition and 
permutation operators on M^. The transposition operator is defined 
by extending the following isometry defined on simple tensors 

U(ij){xi (g) X2 ^ ■ ■ ■ Xi ■ ■ ■ Xj ■ ■ ■ ^ Xn) = Xi® X2® ■ ■ ■ Xj ■ ■ ■ Xi - ■ ■ ® Xn 

to all of Mn. For any vr in the group H^v of permutations of {1,2,..., A^}, 
one may define the permutation operator [7,^ as C7(ij.jj.) • • • ^(i2i2)^{jiji)' where 
{ikjk) ■ • • (*2i2)(*iii) is any product of transpositions that equals vr. 
The antisymmetric subspace may now be defined as 

An = {i^ & IHn : U^iIj = sgn(7r)V' Vtt G Un] ■ 
One may verify that 

is the orthogonal projector whose range is ^at. 

The pure states of an iV-fermion system correspond to the orthogonal 
projectors onto one-dimensional subspaces of An- That is, a pure state 



4 



is given by 

for some iJj G An of unit length. The statistical states of the A^-fermion 
system are the positive trace class operators or density operators D on An 
of trace 1. These can be identified with density operators D on all of JHn 
whose eigenvectors lie in ^at, i.e., such that 

oo 

D = Y^XiP^^ 

i=l 

for some orthonormal system {V'i} in -^A^ and a family of positive numbers 
Aj that sum to 1. It follows that these fermionic densities are those density 
operators that satisfy 

DU^ = U^D = sgn(7r)D Vvr e H^. (2) 

If a density operator D on Mn commutes with every permutation operator 
then it is symmetric. In particular, fermionic densities are symmetric by 

(I)- 

If {ejljgj is an orthonormal basis of JH then 

{eji ® O • • • ejj, : ji, j2, ...Jn^J} 

is an orthonormal basis of Mn- Since is the range of P_Aj^ and since 

-P4iv(eji (gi • • • ejjy) = 

unless all of the indices ji are distinct, the set 

{-P^iv(eji ® (8) • • • Cjj^) : ji, j2, • • • ,jn all distinct} 

is a spanning set for An- In fact it is an orthogonal basis for An, each vector 
having norm 1 / Vectors of the form -v/iVT {ej^ (E) ej^ ^ ■ ■ ■ ^ e^^ ) 
are known as Slater determinants. 

The trace class operators on a Hilbert space M form a Banach space 
T{IH) with the norm ||T||tr = Tr(|T|). The important inequality 

\\TB\\tr < \\T\\tr\\B\\ (3) 

holds whenever i? is a bounded operator of norm ||i?|| and T £ T{]H). It is 
this basic inequality that will produce our key estimates. 



5 



For n < N , the n partial trace is a contraction from T{]H'^ ) onto 
T{H®"-). The n*^ partial trace of T will be denoted Tn, and may be defined 
as follows: Let O be any orthonormal basis of IH. If T € T{IH®^) and 
n < N then 

(T;„(w),x) = 

<^ T(t(; (g) Zi (g) • • • (g) 2:Ar__„), X (g (g • • • (g ZN-ri^ (4) 

2l,...,2jV-nS O 

for any w,x e M^'"-. If a trace class operator T G T{1H^^) satisfies (|) 
then so does r.„, i.e., the partial trace defines a positive contraction from 
T{IH®^) to T[]H®^) that carries fermionic densities to fermionic densities. 

In the following definition, and throughout this article, we use the su- 
perscript ®" to denote the n*^ tensor power of an operator, and we use the 
notation S„ for n!P4^, i.e., 

E„ = ^ sgn(7r)[/^ . 
7ren„ 

The n*^ tensor power of an operator ^ on iff is the operator A®"^ on 
defined on simple tensors by 

A®'^{xi (g X2 • • • (g x„) = Axi AX2 (g • • • (g AXn ■ 



Definition 2.1 For each N , let he a symmetric density operator on 
JHfyf. The sequence {D^} has Slater closure if, for each fixed n, 

lim WDn-u - ^^iSnlL = 0. 

This terminology is motivated by the observation that, if 'I' at is a Slater 
determinant in An and P^^ denotes the orthogonal projector onto the span 
of ^'at then 

7ren„ 

for this implies the following: 

Proposition 2.1 For each N let be a Slater determinant in An, and 
let Pijijy denote the orthoprojector onto the span of ^n- Then {i^jvl ^^■^ 
Slater closure. 
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3 The time-dependent Hartree-Fock equation 



We are going to prove that, in the mean field limit, the time-dependent 
Hartree-Fock equation describes the time-evolution of the single-particle 
state in systems of fermions. We state our theorem in this section and 
go on to prove it in the three subsequent sections. 

First we describe the A^-particlc Hamiltonian. Let L^^^ be a self-adjoint 
operator on M, where L^^^ may depend on N in an arbitrary manner. The 
free motion of the j*^ particle is governed by 

where / denotes the identity operator on JH. The interaction between the 
particles has the form 1/N times the sum over pairs of distinct particles of a 
two-body potential V. Let F be a bounded Hermitian operator on M ® M 
that commutes with the transposition operator C^(i2)- Define the operator 
Vi2 on Mn by 

^12 (xi (g) a;2 (8) • • • <8) xn) = V{xi (g) X2) (8) xa (g) • • • (g) xat 

and for each 1 < i < j < N define Vij = U*Vi2U-,r where tt is any permuta- 
tion with 7r(i) = 1 and 7r(j) = 2. Let 

l<j<N l<i<j<N 

be the N-particle Hamiltonian operator on IHn- The von Neumann equation 
for the iV-particle density operator -Div(i) is 

ihj^D^it) = Yl [Lf\Dr,{t)] + ^ E iV,j,D^{t)] (7) 

l<j<N l<i<j<N 

and has the solution 

^-iH^t/hjj^^QyH^t/h^ (8) 

Next we define the time-dependent Hartree-Fock equation. Let L^^^ and 
V be as above, and let ?7(i2) denote the transposition operator on iff (g) iff. 
The time-dependent Hartree-Fock (TDHF) equation for a density operator 
F{t) on iff is 

ih^^Fit) = [lW, F{t)] + [V, F,-{t)],^^ 

F,-it) = (F(i)®i^(i))(l-C/(i2)) = F(irS2 (9) 
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(the subscript ^ on the last commutator denotes partial contraction). Fol- 
lowing 1^, we define a mild solution of equation @ to be a continuous 
function t i— > F{t) from [0, oo) to the real Banach space of Hermitian trace 
class operators such that 

F{t) = e-i^<"V(0)e+^i^'"' - re-^^'"^[y,F,-(.)].,e+^^^'"'d. 

Jo 

for all t > 0. The results proved in ||5| show that @ has a global mild 
solutionis for each Hermitian trace class initial data F[0). Furthermore, 

F{t) = U*F{0)U (10) 

for some unitary operator depending on t and F{0). In particular, the 
operator norm of F{t) is constant. 

The relationship between the A'^-particle system and the TDHF equation 
is the subject of our main theorem. Recall the definition ( ^.l| ) of Slater 
closure. 

Theorem 3.1 For each N , let D]\[{t) the solution to ^ whose initial value 
Dj\f{0) is a symmetric density. Let F^^\t) be the mild solution of the TDHF 
equation ^ whose initial value is F^'^^O) = D]\[-i{0). 

//{Z)jv(O)} has Slater closure then {Z)jv(t)} has Slater closure and 

hm ||Z)jV:i(t)-F(^^(t)||,, = 

for all t > 0. 

4 Two hierarchies, and their difference 

Consider the A^-particle von Neumann equation (0). From now on we will 
suppose that the initial A^-particle density operator Di\f{0) is symmetric, 
i.e., that 

U:Dn{0)U^ = Dn{0) 

for all vr G Hn- (Recall that, in particular, fermionic densities are sym- 
metric.) The symmetry of the Hamiltonian (^ ensures that D]\[(t) remains 

^The solution obtained in Theorem 4.2 of is indeed defined for all positive times 
because the nonlinearity of TDHF satisfies condition (4.1) of |^ — see Proposition 3.5 
there. 
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symmetric for all t. From (^) and the symmetry of Di^{t), it follows that 
the partial trace DN-n{t) satisfies]^ 

+ E [V^,n+l,DN■.n+l{t)l^ (H) 

l<i<n 



The system of equations (11) for Dj\f:i, Dn-2, ■ ■ ■ , Dn-^n^i together with the 
equation (^) for Djv is called the N-particle hierarchy. For our estimates 
later on, it is convenient to rewrite the equations (jl^) of the hierarchy as 



l<i<n 

+ Sn{t,N,DNm (12) 

with 

en • ) = E [4"^^ • ] 



l<i<j<n 

- ^ E [y^,n+l,DN■.n+l{t)l^ • (13) 
l<j<n 

Next we describe another hierarchy, built from "the bottom up" out of 
solutions to the TDHF equation, in contrast to the hierarchy we have just 
considered, which is built from "the top down" starting with solutions to 
(0). If F is a trace class operator, define = F and 

for n > 1. When F depends on t we write F~{t) instead of F{t)~. The 
notation F^it) has already been used in the TDHF equation (^. 

^In the proof of Theorem we proceed as if Dn was a classical solution of equation 
(|^; the case of general initial data as in the statement of the Theorem follows from a 
density argument and the unitarity of the group generated by the A'^-body Schrodinger 
operator. 
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Proposition 4.1 If F is a classical solution f\ of the TDHF equation 
then 



i=i i=i 



where IZn is defined on trace class operators by TZi{X) = (the zero opera- 
tor) and 



n 

Unix) = X^^+^Y. 



{k,n+l) 



for n > 1. 

Proof: For any trace class operator X, 



(14) 



n 



A:=l 



n n 



St), . 



3=1 k=l 

The first equality in (ITsI ) holds because 



(15) 



k=l 



and the second equality in (15) holds because S„ /^(iH') commutes with 

^Consistently with the previous footnote, we are proceeding as if Fat was a classical 
solution of the case of general initial data as in the statement of the Theorem is 
recovered by a density argument and the continuous dependence on initial data of the 
mild solution to (^: see Theorem 4.1 of [^. 



10 



Yl ^j,n+i- From the TDHF equation (^) we calculate that 



d 



ih[ -F{t) ® 

n 
n 

+ 7^„(F(^)) (16) 



i=i 



n n 

j=i k=i 



(k,n+l) 



• 

:n 



By the identity ([l5|), the last sum in (|T6|) equals ^ [V^>+i, ^n+i(0] ' 

i=i 

proving the proposition. □ 
Now let Df^{t) be a solution of the N-particle von Neumann equation 

(0) and let F{t) be a solution of the TDHF equation (|). For 1 < n < iV 
define the n*'^ difference 

EN,n{t) = DN:n{t)-Fn{t)- (17) 



From the N-particle hierarchy equations ( |12D and ( |13D and Proposition [4.1J , 
it follows that 

d " 

+ f„(t,A^,D^(o)) - 7^„(F(^)) (18) 

for n = l,2, ...,A^ — 1. The characters £■ and TZ were chosen to evoke the 
words "error" and "remainder." Indeed, in the next section we find bounds 
on these error terms under conditions on Di\]{{)) and F{G). The rest of this 
section is devoted to showing how such bounds lead to an upper bound on 
the differences Ep^^t). 

To this end, let us define 

Err(t,A^,n) = f„(t, iV, Z)7v(0)) - 7^„(F(^)). (19) 
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Let C/„ i denote the unitary operator exp X]j=i ) and define 

isometries Z^^^^ on the trace class operators by 

Then ZN,n{t) = ul{l\EN,n{t)) has the same trace norm as EN^nit) and 
satisfies 

, . n 

j^Z^At) = ^7V,n+i(t)]^„-^Z^iyErr(t,iV,n) (20) 

i=i 

for n = 1, 2, . . . , iV - 1. From @ it follows that 

< ||^JV,n(0)||,^ + ^^^^ \\ZN,n+l{s)\l^ds 

+ ^l^\\uif{Evv{s,N,n))l^ds 

for n = 1,2, ... ,7V - 1. Recalling that \\ZN,n+iis)\\tr = \\EN,n+i{s)\\tr and 
that \\uif{Erris,N,n))l^ = ||Err(5,iV, n)\\tr, the preceding inequality be- 
comes 

\\ENAt)\\tr ^ <N,n,t) + 1^ \\EN,n+i{s)\l^ds (21) 

if we define 

1 /■* 

e{N,n,t) = ||^^,n(0)||,^ + j_J^ \\Err{s,N,n)\l^ds. (22) 

Beginning from ( ^T|) and iterating the inequality m times (for some m < 
N — n — 1) we obtain our desired bound on the trace norm of Ei^^t): 




(23) 
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5 Error estimates 



In this section we collect the error estimates that will be used to prove 
Theorem 

If Div(O) is a density operator then the solution D]\f{t) of the N-particle 
von Neumann equation (^) is a density operator for all t > 0, and it is clear 
from (13) that 

\\£n{t,N,DN{Omtr < 3— (24) 

for all t. 

Lemma 5.1 If{Di^} has Slater closure then 



lim WDn: 



0. 



Proof: The trace norm of Dj^.-^ equals the sum of the squares of the 
eigenvalues of D^-i. Since the operator norm of D^-^i equals its largest 

111 /2 

eigenvalue, it follows that ||DAr:i|| < . But Dj^.^ 

whence 



D 



N:l\\tr 



{D^:2-I)fi(/-f/(12))}.^ 



< 



tr 



D 



N:2 



tr 



The Slater closure of {-Dat} implies that the right-hand side of the preceding 
inequality tends to as — > oo. □ 

Lemma 5.2 If F is a density operator then \\F~\\tr < 1 for all n. 

Proof: Since $]„ = (Sn)* = ^(^n)^ commutes with F*^", it follows 
that F- = F^'^T.n = (F®") is a nonnegative operator. Thus, the 
trace norm of F~ equals its trace. This trace is 

jl,..;jneJ 

where {ejjjgj is an orthonormal basis for IH consisting of eigenvectors of 
F. This sum may be taken over distinct indices ji, ■ ■ ■ ,jn £ J, since S„ 
annihilates all tensor products ej^ 



E 



distinct 

ii,...,j„eJ 



■ ej^ with repeating factors, so that 



distinct 

ii,...,j„eJ 
< Tr ( F®") 
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as asserted. □ 
The next lemma provides a bound on the trace norm of the remainder 
term TZn{F) when F is a density operator. The bound is proportional to 
the operator norm of F. 



Lemma 5.3 Let TZn be as in and let F be a density operator. Then 

\\nn{F)\\tr < 2n(n-l)||y|| \\F\\. (25) 
Proof: From (|l^ we see that TZn{F) equals 

I j,k=l ) -n 

Since U(^i^^n+i) commutes with F®"^^^ and since S„ / commutes with 
I]j,fc:jVfc^(fc,n+i)^j,n+i, it follows that 7^„(F) equals 



3^k 



Since the trace norm of a trace class operator equal the trace norm of its 
adjoint, it follows that 



\nn{F)\\tr < 2 



n 

{v,,n+lU^k,n+l){F- d?) F)] 



j,k = l 

j^k 



tr 



(26) 



< 2n{n-l)\\[Vn-l,n+lU^n,n+l){Fn ^F)y 

But one may verify directly that 

{K„i,„+i[/(„,„+i)(F-0F)y^ = (/®"-i®F)K„i,„F-, (27) 
so that, by (^ and Lemma |5.2| , 

n,n+l) 

{F-(S>F)\ < \\F\\ \\V\\ \\F-\tr < \\F\\ \\V\\. 

(. J -.n tr 

Substituting this in iM) yields (O). 



To verify (27), choose an orthonormal basis {ej}j^j for IH and check 
that the operators on both sides of (^) have the same matrix elements 
relative to the basis {ej (8) ej : i, j G j}. □ 
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Let F{t) be a solution of the TDHF equation dS). Since the (operator) 
norm of F{t) is constant, it follows from Lemma 5^ that 



WTZniFmUr < 2n{n-l)\\V\\ ||F(0)|| 

for all t > 0. With this estimate for TZn and estimate (U) for £n, it follows 
that Err(i, N, n) of equation ([l9|) satisfies 

\\ETT{t,N,n)\\tr < 2n2||y|| (^^ + ||F(0)|| 

and e{N, n, t) of equation ( p2|) satisfies 

e{N,n,t) = 2n2||y||^(^l + ||F(0)||) 

+ ||^^.n(0)||,,. (28) 

6 Proof of the theorem 

Equipped with the estimates of the preceding sections, we proceed to the 
proof of Theorem |3.l| . 

So, let us assume that L'Ar(O) is a symmetric density for each and that 
the sequence {Div(O)} has Slater closure. Let Diy{t) be the solution of (^) 
with initial value Dn{0), and let F^^^ (t) be the solution of the TDHF equa- 
tion (|) whose initial value is F^^\o) = Djv:i(0). Let {F^^^}~ (t) denote 
(t)|®"5]„ and let E^n{t) denote the difference between Z)jv.„(t) and 

We have the upper bound ( p3[ ) for the trace norm of EN^nit), into which 
we now substitute the estimates (|28|). In the same stroke, we will replace 
the binomial coefficients ("^^Y^) with the larger quantities (n + k)^/nl and 
we will use the fact that sup |||-£/Af,n+m+i 

(s)||^^} < 2 by Lemma ^ Also, 

let us set T = 2\\V\\t/h. We obtain 



^ m 

■ A:=0 

1 9 

+ + ii^^^^(o)||)t'^+' 

■ A:=0 



+ ^(n + m)"r"^ (29) 
nl 
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for m < N — n — 1. Fix T to be less than 1, i.e., fix t < ^/2||y||. For 
fixed n, consider the hmit of the right-hand-side of (|29| ) as N and m tend to 
infinity. The individual terms (fixed k) tend to 0, for ||f('^''(0)|| tends to 
by Lemma 5J and ||-EAr,n+fc(0)||j,^ tends to thanks to the hypothesis that 
{Dn{0)} has Slater closure (recaU that F^^\o) = DN-.iiO)). On the other 
hand, the series on the right-hand-side of d29) are dominated, uniformly with 



respect to m, by a series that converges absolutely for T < 1. It follows that 



lim = 

if t < /i/2||F||. When n = 1, this shows that lim \\DN:iit) - F^^\t)\ 
and consequently 



tr 



(30) 
= 



lim 

N^oo 



tr 



for n > 1 and t < ?i/2||y||. From ( |30D again it follows that, for any n and 
any t < h/2\\V\\, 



lim 



0, 



tr 



i.e., {Z?Ar(t)} has Slater closure. This proves the theorem up to t = /i/2||y ||. 

Let r = ^/3||y||; the previous argument shows that the theorem holds 
for t e [0, r]. At time r, it is no longer the case that Dn-i{t) = F^^^t). 
However, ||SAr^„+fc(r)||^_^ — > and || ^ as tends to infinity — 

to see this, use (|1§) and the fact that ||F(^)(0)|| ^ 0. An argument nearly 
identical to the one above shows that the theorem holds for t £ [r, 2r] . This 
argument may be repeated to establish the conclusion of the theorem on 
each interval of the form [kr, [k + 1)t] for each nonnegative integer k, and 
hence for all t > 0. □ 



7 Appendix: TDHF equations for wavefunctions 

The time-dependent Hartree-Fock equations were invented by Dirac . He 
wrote the TDHF equations as a system of coupled Schrodinger equations 
for A^ time-dependent electronic orbitals, and thence obtained the TDHF 
equation for the density operator. In this article we have studied the TDHF 
in the latter form: as a von Neumann equation for the density operator. 
This appendix explains how to recast the wavefunction formulation of the 
TDHF equations into the language of density operators used in this paper. 
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The starting point in this discussion is the hnear A^-body Schrodinger 



d ^ 1 

ifi-Q^'^N = -ylZ^-fe^^ + ^v ^ V{xk-xi)^N (31) 

k=l l<k<l<N 

where ^'at = ^'Ar(t, xi, . . . ,xj\f) is the A^-particle wavefunction. (Note that 
the interaction term has been multiphed by 1/{N). This scahng has been 
introduced so that — > oo may yield a mean-field equation for the single- 
particle density, namely, the TDHF equation.) The dynamics defined by 
(HI) is unitary on L'^{{1R'^)^). Therefore, 

|^'Ar(t, xi, . . . , xat)!^ dxi . . . (ixAT = 1 

for all t > if the same equality holds at t = (as is the case if I^'at^ 
is meant to be interpreted as the probability density of the system of 
particles in its configuration space). In the language of (|l|), Lfc = —^Ax^^ 
while Vki denotes the multiplication by V{xk — xi). The TDHF equations 
corresponding to (^) may be written as a system of coupled Schrodinger 
equations for orthonormal orbitals '0i(^) s;), ■02(^5 a:), ... , ipN^t, x): 

ih-Q-'4'k{t,x) = — -Ax'il)k{t,x) +\lJk{t,x)—^ j V{x - z)\'ilji{t,z)\^dz 

1=1 •' 

I ^ r 

/ yix-z)Mt,z)Mt,z)dz (32) 
1=1 •' 

The A^ orbitals remain orthonormal at all times; if 'il)i{t,x),ip2{i-,x), . . . , 



ipN{t^x) is a solution of (32) and 



ipk{^,x)4'i{^,x)dx = 5ki then 



ipk{t,x)iJi{t,x)dx = 6ki for all t > . 
One way to obtain the TDHF equations (p3) from the linear A'^-particle 



Schrodinger equation (31) is to solve a variational problem which would 
lead to ( |3lD if unconstrained, but with the constraint that the A'^-particle 
wave equation remain a Slater determinant at all times |^] . This constraint is 
imposed for the sake of obtaining an computationally amenable approxima- 
tion to (|3l|), and it is not justified on physical grounds. In effect, this paper 
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proves that the constraint maintaining Slater determinants at all times is 
rigorously justified in the mean-field limit. 

To see how the orbital form (|32| ) of the TDHF equations relates to the 
TDHF equation @ discussed in this paper, we shall first rewrite @ as an 
equation for the integral kernel of a time-dependent density operator. To do 
this, we need to know how to translate the partial trace into the language 
of integral operators, for equation (^) involves a partial trace. Let T be a 
trace class operator on L^{1R"^ x iR") having an integral kernel p{x,(^,y,r]) 
with x,y G M"^ and £,,rj G M". The partial trace 

T-m is the operator with integral kernel J p{x, z,y, z)dz. 

We may now convert the TDHF equation (^) into an integro-differential 
equation for a time-dependent integral kernel: Let p = p{t, x, y) be the 
integral kernel of the operator Fit) that appears in (|^. Then F2{t) has 
integral kernel 

p{t, xi,yi)p{t, X2,y2) - p{t, xi,y2)pit, X2,yi) , 
while [y, ^2 has integral kernel 

{V{xi - z) - V{yi - z)){p{t,xi,yi)p{t,z,z) - p{t,xi, z)p{t, z,yi))dz , 



and the TDHF equation (^) in the language of integral kernels is 



d 

ifi-g^p{t,x,y) = - — {/:^x-/\y)p{t,x,y) (33) 



+ j - z) - V{y - z)){p{t, X, y)p{t, z, z) - p{t, x, z)p{t, z, y))dz . 

It may be verified that a solution '>pi(t,x), ip2{t, x), . . . , ipiy{t,x) to the 
orbital form of the TDHF equations (|3^ ) yields a solution p{t, x, y) to the 
integro-differential equation (33) via 

1 ^ 

Pit,x,y) = —^7pkit,x)iJk{t,y)- 

k=l 

The rest of this appendix is meant to serve as a key for reading this 
paper with the Schrodinger wave equation (|3l| ) in mind. 
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To the wavefunction \E'Ar(t, •) solution of the A^-particle Schrodinger equa- 
tion (|3l] ) one associates the operator 0^(1) with integral kernel 



PN{t,Xi, . . .,XN,yi, ...,yN)= 'ifNit,Xi, . . . ,XN)'^N{t,yi, ■ 

The natural Hilbert space M in this context is IH = L^{]R^) 
isomorphic to -L^((iR^)^) through the identification 



■ ,yN)- (34) 

and M^'' 



IS 



JV 



V'l (8) . . . (81 V^AT 



tpk{xk)- 



k=l 



The corresponding representation of the permutation group Hn is given by 
the formula 

(C/^^)(xi,...,X7v) = ^'(3;^-l(i),...,X^-l(Ar)) 

Hence the projection P^^^ is given by 



for ^ G and G L^{{]R^)^ 
the formula 



{PAn'^n)[xi, 



,xn) 



^ X] Sgn(7r)^'Ar(x^(i),...,X^ 

" Trenjv 



(TV); 



A wavefunction ^ is antisymmetric if it is in the image of Pan j equiva- 
lently, if 

^'(x^(i), . . . ,rE^(Ar)) = sgn(7r)^'(3;i, . . . ,XAr) 

for all vr G IItv- If is antisymmetric then the rank 1 orthogonal pro- 
jector P-^j^ with integral kernel as in (34) is fermionic in the sense of (^. 
Property @) extends to all convex combinations of such projectors Pm^^, the 
fermionic density operators discussed in this article. Property (^) implies 
that fermionic density operators commute with all of the operators C/tt, so 
that the integral kernel of a fermionic density operator is symmetric in the 
sense that 

/)Ar(xi,...,XAr,yi,...,yAr) = Pn{x^(i), ■ ■ ■ ,X^(N),y-K{l), ■ ■ ■ ^y-K{N))) 
for all vr G IItv- 

In the case where ^ = tpi® . . .®ipjs[ with , . . . , ^/^tv orthonormal, one 
finds that the Slater determinant \/]vT Pan{^n) truly is a determinant: 



N\{PA^'i>N){xu...,XN) 



■0l(xi) i'l{x2) 
^2{xi) 1p2{x2) 

V'iv(xi) Viv(x2) 



■^i{xn) 

^2{xn) 
tpN{xN) 



19 



Acknowledgement . 

This research was supported by the French- Austrian "Amadeus" program 
(OAD V6) and by the Austrian START project "Nonlinear Schrddinger and 
quantum Boltzmann equations" of N.J. M. Also, F.G. was supported by the 
Institut Universitaire de Prance. The first three authors thank the ESI in 
Vienna for its hospitality. We also express our gratitude to Xavier Blanc, 
Eric Cances and Claude Le Bris for numerous and helpful discussions on 
this paper. 



20 



References 



[I] R. Alicki and J. Messer. Nonlinear quantum dynamical semigroups for 
many-body open systems. Journal of Statistical Physics 32 (2): 299 - 
312, 1983. 

[2] C.Bardos, F. Golse and N.J. Mauser. Weak coupling limit of the N- 
particle Schrodinger equation. Mathematical Analysis and Applications 
7 (2): 275 - 293, 2000. 

[3] C. Bardos, L. Erdos, F. Golse, N.J. Mauser and H.-T. Yau. Deriva- 
tion of the Schrodingcr-Poisson equation from the quantum A^-particle 
Coulomb problem. C. R. Acad. Sci. Ser. I Math 334: 515-520, (2002). 

[4] L. Boltzmann. Lectures on Gas Theory. Dover Publications, New York, 
1995. 

[5] A. Bove, G. Da Prato, and G. Fano. An existence proof for the 
Hartree-Fock time-dependent problem with bounded two-body inter- 
action. Communications in Mathematical Physics 37: 183-191, 1974. 

[6] E. Cances, and C. Le Bris. On the time-dependent Hartree-Fock equa- 
tions coupled with a classical nuclear dynamics. Math. Models Methods 
Appl. Sci. 9 (7): 963-990, 1999. 

[7] P. Dirac, Note on exchange phenomena in the Thomas atom, Proc. 
Cambridge Philosophical Society 26 (1930) pp. 376 - 385 

[8] N. G. Duffield and R. F. Werner. Local Dynamics of Mean-Field Quan- 
tum Systems. Helvetia Physica Acta 65, 1016-1054, 1992 

[9] L. Erdos and H.-T. Yau. Derivation of the nonlinear Schrodinger equa- 
tion with Coulomb potential. Preprint (2002). 

[10] A. D. Gottlieb. Markov transitions and the propagation of chaos (PhD 
Thesis). Lawrence Berkeley National Laboratory Report, LBNL-42839, 
1998. 

[II] F. A. Griinbaum. Propagation of chaos for the Boltzmann equation. 
Archive for Rational Mechanics and Analysis 42: 323-345, 1971. 

[12] M. Kac. Foundations of kinetic theory. Proceedings of the Third Berke- 
ley Symposium on Mathematical Statistics and Probability, Vol III. Uni- 
versity of California Press, Berkeley, California, 1956. 



21 



[13] M. Kac. Probability and Related Topics in Physical Sciences. American 
Mathematical Society, Providence, Rhode Island, 1976. 

[14] H. P. McKean, Jr. A class of Markov processes associated with nonlinear 
parabolic equations. Proceedings of the National Academy of Science 56: 
1907-1911, 1966. 

[15] S. Meleard. Asymptotic behavior of some interacting particle systems; 
McKean- Vlasov and Boltzmann models. Lecture Notes in Mathematics, 
1627. Springer, Berlin, 1995. 

[16] H. Spohn, Kinetic Equations from Hamiltonian Dynamics. Rev. Mod. 
Phys. 53: 600 - 640, 1980. (Sec Theorem 5.7) 

[17] A. Sznitman. Equations de type dc Boltzmann, spatialement ho- 
mogenes. Zeitschrift fiir Wahrscheinlichkeitstheorie und verwandte Ge- 
biete 66: 559-592, 1984. 

[18] A. Sznitman. Topics in propagation of chaos. Lecture Notes in Mathe- 
matics, 1464. Springer, Berlin, 1991. 



22 



Resume 



On montre dans ce travail que les equations d'evolution de Hartree-Fock 
decrivent la limite de I'equation de Schrodinger a N corps pour N tendant 
vers I'infini et une constante de couplage en 0{1/N) et pour des donnees 
initiales prochcs de determinants de Slater. On nc considere ici que le cas de 
potentiels d'interaction binaires, symetriques et bornes. Lorsque oo, 
on montre que la suite des traces partielles "a un corps" de I'operateur 
densite a N corps converge, au sens des operateurs a trace, vers la solution 
de I'equation de Hartree-Fock sous forme operatorielle. 

Mots-cles: Equations de Hartree-Fock; Probleme a N corps quantique; 
Approximation de champ moyen 
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